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Abstract 

In a bounce-free spherical hydrodynamic implosion, the post-stagnation hot core plasma does 
not expand against the imploding flow. Such an implosion scheme has the advantage of improving 
the dwell time of the burning fuel, resulting in a higher fusion burn-up fraction. The existence of 
bounce-free spherical implosions is demonstrated by explicitly constructing a family of self-similar 
solutions to the spherically symmetric ideal hydrodynamic equations. When applied to a specific 
example of plasma liner driven magneto-inertial fusion, the bounce-free solution is found to produce 
at least a factor of four improvement in dwell time and fusion energy gain. 



PACS numbers: 52.58.Qv, 52.35.Tc, 52.50.Lp, 52.25.Xz 



Spherical hydrodynamic implosion is a central concept in inertial confinement fusion 
(1CF) pQ and magneto-inertial fusion (MIF) [2J, which use a variety of inertial pushers 
driven by lasers [3], heavy ion beams [J], and pulsed power It is also of fundamental 
importance in late-stage stellar evolution, especially supernova core collapse [6]. In a typical 
spherical hydrodynamic implosion, converging mass flow provides compressional (and pos- 
sibly shock) heating of the core, converting flow kinetic energy into thermal energy of the 
plasma localized about the center of symmetry. Once the peak pressure is reached, a stag- 
nation surface appears between the core and imploding flow plasmas. What happens next 
is usually a bounce motion of the core in which the hotter plasma at the center expands 
against the imploding flow. In the special case where the interface between the hot core 
and imploding flow remains stationary rather than expanding outward, one has a bounce- 
free spherical hydrodynamic implosion. That is, while in the typical bounce case both the 
core/pusher interface and stagnation surface move radially outward, in the bounce-free case 
the core/pusher interface remains stationary and only the stagnation surface moves radially 
outward. In either case, once the stagnation surface sweeps through the entire imploding 
material, the whole core-pusher assembly will undergo an expansion. The latter imposes an 
upper limit on the system dwell time. 

Optimizing toward a bounce-free implosion provides the most design freedom in increasing 
the dwell time and thereby a higher fusion burn-up fraction, which as shown in this Letter, 
can be achieved by employing a specially shaped pusher. This approach is most naturally 
applied and would be of most benefit to plasma liner driven MIF [7] since there the pusher 
(liner) is created by merging plasma jets conceptually capable of modulating the imploding 
flow profile. Also, due to the challenge for MIF to attain an ignited burn wave, it has the 
greatest need for increasing the burn fraction via increased dwell time. Conversely, the lack 
of a strong burn wave and intense hot spot pressures might actually make the bounce-free 
implosion scheme energetically feasible to implement. 

In this Letter we present a family of self-similar solutions to the spherically symmetric 
ideal hydrodynamic equations which produce bounce-free spherical implosions. The physics 
implications of these bounce-free implosion solutions are explained in the application to 
plasma liner driven MIF. In particular, we find that by maintaining a bounce-free implosion 
through liner (inertial pusher) profile shaping, the dwell time of the compressed hot fuel can 
be improved by at least a factor of four compared with a specific example of a stationary 
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liner profile [8], while holding the liner speed and total energy the same. 

To motivate the analytical solution, we first outline the physical picture of and the math- 
ematical constraints on bounce-free spherical hydrodynamic implosions. Figure [T] shows a 
sketch of the assembly after the liner hits the target. The stagnation or ignition pressure of 
the target p s t is achieved at time t st with the spherical interface between the hot spot and 
the liner at r st coming to rest. To be bounce-free, this interface remains still for t > t st . 
This is made possible by shock heating of the liner plasma neighboring r st . In other words, 
a shock arises at the interface r st at t st and propagates radially outward. The shocked liner 
in the region of r st < r < R s ^ is in static pressure balance with the target, p(r, t) = p st . The 
shock front at i? s h(£) becomes a second stagnation surface, which along with the interface 
at r st , bounds an ever-expanding, stagnated, and shock-heated liner. The imploding liner 
for such a bounce-free regime must have highly constrained profiles of density p, pressure p, 
and flow speed u. 

This Letter demonstrates the existence of bounce-free spherical implosions by presenting 
a family of imploding flow solutions that exactly satisfy the hydrodynamic equations and 
jump conditions at R s h(t) as well as the following constraint 

p(r,t) = const., p(r,t) = const., and u(r,t) = (1) 

for < r < R s h(t) and t > t st . Interestingly, such a solution consistently describes cases 
of both an infinitesimally small target (void) and a finite sized target kept still. Indeed, 
the liner flow solution with no distinct target satisfying Eq. ([I]) obeys all the first principle 
equations even upon "inserting" the target into the center of symmetry, while condition of 
the shocked liner pressure being constant ensures that the target stays at rest. We notice 
that the solution we are aiming for is thus a generalization of a well known Noh's solution [H] . 
Indeed, Noh's solution gives the uniform pressure and density in the stagnated gas for the 
case, in which the initial converging shell, and therefore the the unshocked part of the liner 
as sketched in Figure [TJ is perfectly cold, while the solution to be presented here allows 
temperature of the shell to be finite. Unlike Noh's case, it is then not possible to write a 
simple analytical formula describing liner evolution. It is still possible though to reduce the 
problem to a set of relatively simple ordinary differential equations, as we demonstrate later 
in this Letter. 

The simplest model for spherically symmetric ideal hydrodynamic implosion is the con- 
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FIG. 1: Illustration of a bounce- free spherical implosion. The target (in gold) stays still within 
its boundary located at the stagnation radius r s t- The shocked part of the liner (in red) is at rest 
with pressure and density being uniform and equal to p s t and p s t respectively. The shocked front 
(denoted by R s h) is propagating outward, while the unshocked liner (in blue) is streaming inward. 

servation laws for fluid mass, momentum, and entropy 

(pur 2 ) 

Pt + ^^ = 0, (2) 

U t + UU r + — = 0, (3) 
P 

{pp-^) t + u(pp-^ r = 0, (4) 

where the subscripts "t" and V denote differentiating over the time and radial variables, 
respectively. These three equations need to be solved in the unshocked region of the liner 
with the jump condition at the shock front R s h{t) connecting the infalling liner to the 
bounce-free constraints, Eq. Q. To this end, it is convenient to focus our consideration on 
the similarity solutions. A large class of such solutions and their application to ICF problems 
is considered in Ref. [3]. In our case, its sub-class, the so-called spherical quasi-simple waves, 
turns out to be sufficient. To relate these waves to those described in Ref. |3], one should 
start from the more general form given by its Eqs. (6.182-6.183) and set a = 1 and k = 0. 
In the next paragraphs we follow §§162-164 of Ref. pU] to outline the basic features of this 
solution family which we then utilize for the purpose of this Letter. 

Spherical flow belongs to the class of quasi-simple waves if the velocity, pressure, and 
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density are constant on rays rj = t/r — const. An immediate consequence of this property 
is that the flow is isentropic, i.e., pp -7 is not only conserved for a given fluid element, as 
guaranteed by Eq. Q, but is also the same for all fluid elements. The pressure and density 
can be then expressed in terms of the sound speed c and with the help of some algebraic 
manipulations, Eq. ^ and (|3| reduce to 

DrjU v = \(1 - U) 2 - 3C 2 ]u, (5) 
Dr)C v = [(1 - Uf - (7 - 1)17(1 -U)- C 2 ] C, (6) 

where D = (1 — U) 2 — C 2 , U = 7]u(t]) and C = rjcij]). To write Eqs. (5]6), pressure, 



density, velocity, the radial coordinate and time are normalized to p s t,Pst, ^Pst/Pst, r st and 



r s t/\/p s t/p s t, respectively. 

It is instructive to consider possible scenarios on the (U, C) plane. To do so we divide 
Eq. (§) by Eq. Q to find 

dC = (l-y) 2 -( 7 -l)y(l-y)-C 2 C 

dU (1 -U) 2 -3C 2 V 1 ' 

The full vector field of Eq. ^ is quite complex and contains a number of different regimes 
(e.g., see Fig. 11 of Ch. 6 of Ref. [10J ) . The C{U) trajectories of interest (i.e., those cor- 
responding to the converging flow stopped by the shock wave propagating outward) are 
located in the upper-left quadrant of the plane, as plotted in Fig. [2] The arrows are in the 
direction of growing 7], i.e., the direction of growing t for fixed r or decreasing r for fixed t. 

The jump conditions along with Eq. can be employed to find that C s h and U s h, values 
of C and U just before the shock front, must obey 

(1 - fi 2 )C 2 h = (1 - p 2 )(l - U sh ) 2 + U sh (l - U sh ) (8) 

c/ sh = (i-/i 2 )(i-4K), (9) 



where t> s h is the shock speed and /i 2 = (7 — l)/(7 + 1). Eq. (J8J) implies that the shock 
transition takes place on an ellipse in the (U, C) plane. Choosing a certain point on this 
ellipse then gives the shock speed u s h through Eq. Hence, the desired regime is described 
by a curve uniquely defined by f s h that starts somewhere on the ellipse in the upper-left 
quadrant and runs toward the origin of the (U, C) plane, since this direction corresponds 
to going to larger radii for a fixed t. Importantly, C and U tend to zero together with 77, 
so that the sound speed c = ri~ l C{rj) and the fluid velocity u = rj^Uirj) stay finite. For a 



given curve, U and C as functions of 77 can be found by integrating Eqs. ( 8~p ) from 77 = 77 s h 
to rj = with initial conditions U(rj sh ) = f/ sh and C(r] sh ) = C s h, where r] sh = l/v sh . Once 
U{rf) and C{rj) are evaluated, w(?y) and c(r]) can be recovered from definitions given after 
Eq. (|6j), and p(^) and 79(77 ) can be found from the isentropicity condition. As a result, we 
obtain a solution family parameterized by v s h for bounce-free spherical imploding flows. 

Any set of profiles v Vsh (rj) and c Vsh (rj) generated as described in the preceding paragraph 
explicitly answers the question of what the liner profiles must be at and after stagnation 
in order for the target to maintain constant pressure and remain bounce-free. Next, we 
investigate the evolution of the liner prior to stagnation to see whether the liner profile 
found can be created in a practical device. To do so, we need to solve Eq. ( 6p ) backward 
in time, i.e., extend the previously found solution from the upper-left quadrant of Fig. [2] to 
the lower-right one. In physical space, passing through the origin of Fig. [2] corresponds to 
the liner reaching the center of symmetry. Consequently, unlike the previously considered 
post-stagnation solutions, this part of our analysis (for t < 0) only qualitatively applies to 
the physical case of a finite-size, preformed target. Nevertheless, even such a qualitative 
consideration provides substantial insight into relevant liner parameters, as we demonstrate 
next. 

As shown in Fig. [2j a critical point exists on the ellipse such that for v S h > v CT the solution 
curve ends up in the "rest point" (U — 0, C — — 1), whereas for v 8 t < v CT in the "cavitation 
point" (U = 1, C = 0). The former case corresponds to the situation in which initially the 
entire space beyond the stagnation radius is occupied by a liner material at rest with finite 
density and pressure. The latter case corresponds to the situation in which initially there is 
a vacuum between the infalling liner material and the target. This is the default regime for 
a plasma liner driven MIF experiment. Despite behavior of all the curves being qualitatively 
the same for 77 > 0, only those generated from < v CT shall therefore be considered. To 
show the characteristic liner evolution in the cavitation regime we set t> s h = 0.88 < v CI and 
integrate the fluid equations as previously described to obtain v Vah (r]) , c Vah (rj) and p Vah {v)- 
Then, by inserting t/r for 77 we generate a series of snapshots to create an animation [TT] of 
pre- and post-stagnation phases of a bounce-free spherical implosion. 

Having now proven the existence of bounce-free spherical hydrodynamic implosion by 
explicitly constructing a family of self-similar solutions to the spherically symmetric hydro- 
dynamic equations, our next objective is to elucidate the physics implications of this regime. 
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FIG. 2: Integral curves of Eq. for 7 = 5/3. The upper-left quadrant (77 > 0) corresponds to 
the bounce-free implosion regime sketched in Fig. [TJ The red branches in the lower-right quadrant 
(?7 < 0) correspond to liner convergence in vacuum before reaching the origin. Trajectory in green 
corresponds to Noh's solution. 

As previously mentioned, the physics advantages of bounce-free implosion is probably most 
easily appreciated in the plasma liner driven MIF concept. This is for two reasons. The first 
is that generally MIF schemes do not rely on an ignited burn wave because they operate at 
ignition densities substantially lower than for ICF; the desired high burn fraction is primar- 
ily achieved by maintaining the compressed hot spot at thermonuclear burning condition 
for as long as possible. A bounce-free implosion minimizes the decompressional cooling of 
the target compared with a normal bounced implosion. The second is that by employing 
a standoff driver such as an array of plasma guns to form and deliver the liner, one has 
substantial freedom in shaping the liner profile (density, speed, and thermal pressure) and 
the possibility of using a thick, as opposed to thin, liner. It is also of interest to note that 
because of the much smaller areal density of MIF targets, compared with that of conven- 
tional ICF, the mean free path of an alpha particle is much greater than the target radius 
and neglecting fusion heating effect on the target pressure and on the right side of Eq. Q 
is a useful approximation. Next, we illustrate the potential energy gain by comparing the 
dwell time estimate of the bounce-free implosion solution with an example of a conventional 
implosion. 

Discussion on which of the curves with v s h < v CI would work best for a realistic plasma 
liner MIF device is beyond the scope of this Letter and will be presented elsewhere. In 
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what follows, we focus on the general effect of liner shaping toward bounce-free implosion 
to demonstrate that employing a solution presented here can result in a substantial dwell 
time improvement. To do so, it is convenient to take the estimate of Ref. [5], where no 
special liner shaping is assumed and a simple stationary flow solution is used to model the 
liner profiles, as a reference point. Namely, we next consider a liner of the newly found form 
whose energy and velocity, as well as the relevant target parameters, are the same as in 
Ref. j8] . It is then reasonable to attribute the difference between the dwell time predictions 
to that between the liner shapes and, accordingly, whether or not the implosion results in 
target bounce. 

The earlier estimate [Hj assumes that a target expands (bounces) at a speed as fast as the 
shock front propagation for t > t st , and the fuel disassembly is considered to be complete 
once its radius becomes 2r st as given in Eq. (38) of Ref. [8]. The dwell time is therefore 
obtained by computing r st /t> sh , where r st is 0.5 cm. The shock speed is estimated to be half 
of the 100 km/s liner velocity and is thus 50 km/s. As a result, the dwell time estimate is 
100 ns. The target temperature is taken to be 10 keV, with a fuel burn-up fraction of about 
0.01, and the corresponding stagnation pressure is 62.5 Mb. The total liner kinetic energy is 
122 MJ and the fusion energy gain factor G is estimated to be about 2.6%. We now proceed 
by conducting a similar estimate for the bounce-free liner solution. 

If it were possible to maintain the unshocked liner flow with the prescribed profiles in- 
finitely long, the target would never bounce and be forever kept in a state of maximum 
compression. Practically, the total liner energy is finite and we need to "cut" the liner pro- 
file at a certain point to match the above-mentioned 122 MJ total liner energy constraint. 
Of course, such cutting of a fluid equation solution results in a rarefaction wave spreading 
beyond the outer boundary of the liner. However, until a moment £*, when the shock wave 
hits the outer boundary of the liner, this rarefaction wave propagates at the unshocked liner 
sound speed, which is generally much less than the speed of the liner itself. It is therefore a 
good approximation that despite being finite the liner evolves in space and time according to 
the exact solution above for t < £*. At t > t*, the entire liner is shocked and the rarefaction 
wave is no longer negligible, since the shocked liner sound speed c st is comparable to v^. 
Once the rarefaction wave hits the target, it finally starts disassembling. Accordingly, the 
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dwell time Td w can be estimated by 



r dw — (t* ~ t s t) + ^rare; (10) 

where t rare is the time it takes the rarefaction wave to travel from R*, the outer liner boundary 
position at t = t*, to r st . Note that we are neglecting the further finite time it would take for 
the target to double its radius, which is the criterion for dwell time used in Ref. [S]. Thus, 



our estimate, Eq. (10), is a conservative one. 



To evaluate the first term on the right hand side of Eq. ( 10 ) we note that t* — t st = (R* — 
r s t)/v s h, where i?* can be calculated from liner energy conservation since the hydrodynamic 
efficiency of the target compression is low [HJ [12] . That is, 

MR* - rj) Pst _ p ( , u 

o 7 — Airier, \ Li -> 

o 7 — 1 

where we have used the fact that the energy density of the shocked liner is equal to p s t/(7~ 1) 



and -Euner stands for the total liner energy. It should be noted that Eq. (11) gives R* oc 



(7 — 1) 1//3 , making our estimate quite insensitive to uncertainties in 7. To isolate the effect 



of bounce-free liner shaping, we insert into Eq. (11) the same values of 62.5 Mb, 0.5 cm, 
and 122 MJ for p s t,r st , and -Er me r as in Ref. [SJ. Eq. (11) then gives R* ~ 1.5cm = 3r st and, 
upon setting the liner and shock front velocities 100 km/s and 50 km/s to match those of 
Ref. [8], we find U - t st « 200 ns. 



The second term on the right hand side of Eq. ( 10 ) can be estimated by noticing that 
the rarefaction wave propagates at c st ~ f s h, where, as in the preceding paragraph, we set 
f s h = 50 km/s. It then takes about (R* — r st )/w S h ~ 200 ns for the rarefaction wave to reach 



the target and Eq. (10) gives r^ w ~ 400 ns, which is four times larger than its counterpart 
in the bounced implosion. Since the target parameters are chosen to be the same, the fusion 
energy gain factor is four times larger as well (even neglecting the target radius doubling 
time in the bounce-free case). 

In conclusion, we have demonstrated that by shaping the profiles of an imploding iner- 
tial pusher, the concept of bounce-free spherical hydrodynamic implosion can be physically 
realized. A family of self-similar bounce-free solutions to the spherically symmetric ideal 
hydrodynamic equations is explicitly found, along with a description of their experimental 
accessibility. For the specific application of plasma liner driven MIF, we show that using 
bounce-free liner profiles can substantially slow down the core plasma expansion. As a result, 



the dwell time, and therefore the fusion gain, can be noticeably increased over an unshaped 
liner with equivalent kinetic energy. Compared to the specific case treated in Ref. jH], em- 
ploying the bounce-free implosion regime improves the dwell time and energy gain by at 
least a factor of four. Furthermore, the newly found implosion regime supports the idea 
of using deuterium-tritium fuel in the inner parts of the liner (as suggested in Ref. [7j and 
Ref. [2]), which upon becoming shocked will also burn, thus further increasing the gain. 
Indeed, in such a regime, the post-shocked liner is at rest, i.e., in contrast to previously 
considered schemes, the kinetic energy of the original liner is entirely converted into internal 
energy. This feature brings the temperature of the liner next to the target closer to fusion 
relevant magnitudes, which may further improve the overall efficiency of plasma liner driven 
MIF devices. 
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